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ABSTRACT 

A model for a channel is given. For this model, the 
recursive method is presented in order to calculate the 
probability of K symbol errors in a block of n m-bit symbols. 
The blocks can be interleaved or not. 
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SECTION I 
THE MODEL 

In estimating the performance of an error-correcting device on 
a specific communication channel, it is necessary to find a meaning- 
ful, yet tractable, mathematical model for that channel. An examination 
of data from real channels suggests that most channels pass through 
three distinct phases. The first phase, which nearly any error- 
correcting scheme can handle successfully, is that of long periods of 
practically error- free transmission. The second phase, which is the 
antithesis of the first phase in that no scheme can expect to correct 
it , is that of complete loss of signal for substantial periods of time. 
The third phase might be described as a generalized sputter or bursts 
of error bursts, and it is this phase, if it occurs frequently, for 
which error correctors should be designed. This third phase may 
sometimes be described by means of a two- state model. The three 
phases suggest a Markov process with four states, but such a process 
is mathematically unwiedly. However, when the three phase picture is 
reasonably correct and the sputter phase is accurately modeled for 
purposes of estimating coder performance, it is a simple matter to 
make corrections in such estimates for those periods when transmission 
is nearly error free or when the signal is lost. 

The use of such a two- state model was suggested to us by some work 
by Gilbert [ 1] . Gilbert describes a model for binary error distributions 



in channels subject to noise bursts. Let {x.} be the error process with 
x. = 1 for an error in the i demodulated bit, x. = for no error. 

Two states, designated G and B, of the channel are postulated such that 

t- v» 
at the i bit the state S. is G or B, and the state S. ( , at bit: i + 1 

depends only on S.. Thus, the state sequence {S.} is a simple Markov 

chain described by the two transition probabilities, P: G -» B and 

p: B -* G. We also use Gilbert's notation Q = 1 - P and q = 1 - p for 

the G -> G and B -> B transitions respectively, Let h,k denote respectively 

the probabilities of a correct bit in B and in G Then P{x.=l} = 1 - h = h ' 

if S. = B and Pfx.-l] = 1 - k = k 1 if S. - G. 
l l l 

Using this model with k = 1, Gilbert obtained a good fit to certain 
phone line error data. The statistic fitted was the probability of 
occurrence of zero (= error free) runs of length at least K (The use 
of k < 1 has also been considered by Elliott [5].) 

At MITRE, we have found sets of values for the model parameters 
P, p, h, k, which yield good fits to error data for several different 
types of communication media. One important statistic involved is the 
probability of specific error densities in various length blocks. 

Given that we have p, P, h,k, [2,3,4], we present the recursive 
method used to calculate the probability of K symbol errors in a block 
of n symbols where each symbol consists of m bits. 



SECTION II 
RECURSIONS - PART I 

In Appendix II of [2] (or Appendix B of [3] or [4]) we presented a 
brief outline of the recursive method used to calculate the probability 
of K symbol errors in a block of n symbols where each symbol consists 

of m bits. Now we present the outline in full. Since the last 
documentation, the technique has been extended to permit the n-symbol 
blocks to be interleaved or time spread. (To interleave s blocks means 
to transmit sequentially the first symbol of each of s blocks followed by 
the second symbol of each of those blocks, etc. Thus, a given m-bit 
symbol is transmitted s symbol times after the symbol preceding it in 
its block. ) 

Letting T and U represent either of the states G and B, we define 

TOU(t) = P(x = • • • = x * and S - U|S - T) 

TlU(t) - P(for some i £ t, x =1 and S * U|S Q = T) 

Then 

G0G(1) = Qk G0B(1) = Ph 

G1G(1) - Qk 1 G1B(1) - Fh 1 

BOB(I) « qh B0G(1) « pk 

B1B(1) - qh f B1G(1) - pk f 

and 

GOG(t) = [GOB(t-l) p + GOG(t-l) Q] k 

GlG(t) * [GOB(t-l) p + GOG(t-l) Q] k f + GlB(t-l) p + GlG(t-l) Q 

3 



GOB(C) = [GOB(t-l) q + GOG(t-l) P] h 

GlB(t) = [GOB(t-l) q + GOG(t-l) P] h' + GlB(t-l) q + GlG(t-l) P 

BOB(t) = [BOB(t-l) q + BOG(t-l) P] h 

BlB(t) = [BOB(t-l) q + BOG(t-l) P] h' + BlB(t-l) q + BlG(t-l) P 

BOG(t) = [BOG(t-l) Q + BOB(t-l) p] k 

B1G(C) = [BOG(t-l) Q + BOB(t-l) p] k 1 + BlB(t-l) p + BlG(t-l) Q 

Shortly after the program was written, we discovered that GOG(m), 
GOB(m), GlG(m), GlB(m), BOG(m), BOB(m), BlG(m), BlB(m) could be obtained 
from a difference equation in powers of J and L (see below). Since on 
our computer (IBM 7030), it took under a second to compute all eight 
quantities, we decided not to use the difference equation. However, we 
work out two, and give all eight results. 

GOG(t) = [GOG(t-l) Q + GOB(t-l) p} k 
GOB(t) = [GOB(t-l) q +G0G(t-l) P} h 

The eigenvalues come from the 2nd order linear difference equation: 

f t+l " (Qk + qh) f t " (P_Q) f t-l = °' 
that is, 2J = Qk + qh + V(Qk + qh) 2 + 4hk(p-Q) 

and 21, = Qk + qh - \/(Qk + qh) 2 + 4hk(p-Q) . 
Thus we have 



GOG(t) = Q^ J t + o^ L* 
GOB(t) = /3 L ^ + /3 2 h 1 



and we get q^ , q^ , ft and fl 

from the initial conditions 

GOG(O) = 1, G0G(1) = Qk 
and GOB(O) = 0, G0B(1) = Ph . 

So of, + OU =" 1, Of, J + OUl = Qk, which yields 

GOG(t) = {(Qk - L)/(J - L)) / + {(J - Qk)/(J - L)) L fc 

Also, + j3 2 = 0, ^J + 2 L = Ph, which yields 

GOB(t) = {Ph/(J - L)} (/ - L fc ) 
All eight solutions are as follows: 

GOG(m) = {(Qk - L)/(J - L)} J™ + { (J - Qk)/(J - L)) L™ 

GOB(m) = {Ph/(J - L)) (J™ - L m ) 

BOG(m) = [pk/(J - L)} (J™ - L m ) 

BOB(m) = {(qh - L)/(J - L)) J™ + {(J - qh)/(J - L)) L™ 

GlG(m) = p/(p + P) + [P/(p + P)) (Q-P) m - {(J - qh)/(J - L)) J™ 

- {(qh - L)/(J - L)) L m 

GlB(m) = — [1 - (Q-p) ] - — (J - L ) 



BlG(m) - -^ [1 - (Q-P)'"] - *TI ( J " L ) 



SECTION III 
RECURSIONS - PART 2 

Again letting T and U represent either of the states G and B, 
we define 

TOUI(s) = P(x, , . .. = x, - . = ■ - • = x = 0,S = u|S n = T) 
(s-l)mfl (s-l)mf2 sm sm '0 

* P(m-bit symbol after s symbol times is correct and 

ends in state UJstate T) 

TlUI(s) - P(for some 1 £ i £ m, x, . . , . = 1, S =u|S^=T) 

(s-l)mfi sm '0 

* P(m-bit symbol after s symbol times has at least one bit 

error and ends in state u|state T) 

Let GXG = GOG(m) + GlG(m), GXB - GOB(m) + GlB(m), BXG « BOG(m) + BlG(m), 
and BXB = BOB(m) + BlB(m). 

Then GOGI(s) = GXG • GOGI(s-l) + GXB ' BOGI(s-l) and BOGI(s) = 
BXG GOGI(s-l) + BXB • BOGI(s-l). 

(There will be similar equations in G1BI, B1BI and GOBI, BOBI , and 
G1BI, B1BI, but they have the same eigenvalues and will be omitted.) 

TOGI(s) - [GXG + BXB] TOGI(s-l) + [GXG-BXB-GXB'BXG] T0GI(s-2) - 

This yields the eigenvalues: 



2 a = GXG + BXB + ylGXG + BXB] 2 -4 [GXG • BXB - GXB • BXG] 



2 t = GXG + BXB - ylGXG + BXB] -4 iGXG ' BXB - GXB • BXG] 



Since GXG + GXB = BXB + BXG = 1, 

<j"t = GXG • BXB - GXB • BXG = (1-GXB) (1-BXG) - GXB • BXG = 1-GXB-BXG 

a + t = GXG + BXB = 2 - GXB - BXG = 1 + a T . 

Hence a « 1, and T = 1 - GXB - BXG = GXG + BXB - 1. 
Thus we have 

GOGI(s) = \ • 1 S + X 2 T S 

and BOGI(s) = u * I s + /i T S 

and we get X , X , fi and ji from the initial conditions. 

G0GI(1) = GOG(m), G0GI(2) = GXG - GOG(m) + GXB • BOG(m) 
BOGI(l) = BOG(m), B0GI(2) = BXG * GOG(m) + BXB • BOG(m). 

Hence 

X ' 1 + X ' T - GOG(m), X • l 2 + X £ . • T 2 = GXG • GOG(m) + GXB • BOG(m). 

Solving, X = [GOG(m) (1 - BXB) + BOG(m) • GXB] / (1 - t) 

\ = [GOG(m) BXB - BOG(m) GXG'] I (1 - t) - [GOG(m) - X ] / T - 

Also, 

11 • 1 + ^ • t - BOG(m), ^ • l 2 + fy ' T 2 = BXG * GOG(m) + BXB - BOG(m). 

Solving, ii - [BOG(m) (1 - GXG) + BXG * GOG(m)] / (1 - t) 
^ - [BOG(m) - M x ] / T 



ven 



Therefore, GOGI(s) - X * 1 + X 2 • t , and BOGI(s) - /x • I s + j^'T , 
where X- , X~ , jll , /!« are given above. Since s is fixed for any gi 
application, we will refer to GOGI(s) as GOGI, and to BOGl(s) as BOGI. 
Consider the first i m-bit symbols of a random interleaved block. 
Let GB(i,j) = P(j symbol errors in i symbols and S. * B |s n * G). 
Similarly, we define GG(i,j), BB(i,j) and BG(i,j). 
Then 

GG(1,0) « GOGI GB(1,0) - GOBI 

GG(1,1) - G1GI GB(1,1) = G1BI 

BG(1,0) = BOGI BB(1,0) * BOBI 

*BG(1,1) = B1GI BB(1,1) = BlBI 

Finally, for i = 2 , • • ■ , n and j =0, ' " , i 

GG(i,j) = GG(i-l,j) GG(1,0) + GB(i-l,j) BG(1,0) 

+ GG(i-lJ-l) GG(1,1) + GB(i-l,j-l) BG(1,1) 

GB(i,j) = GG(i-l,j) GB(1,0) + GB(i-lJ) BB(1,0) 

+ GG(i-l,j-l) GB(1,1) + GB(i-lJ-l) BB(1,1) 

BG(i,j) - BG(i-lJ) GG(1,0) + BB(i-lJ) BG(1,0) 

+ BG(i-lJ-l) GG(1,1) + BB(i-l,j-l) BG(1,1) 

BB(i,j) = BG(i-l,j) GB(1,0) + BB(i-l,j) BB(1,0) 

+ BG(i-l,j-l) GB(1,1) + BB(i-l,j) BB(1,1). 



Finally, 



P(random bit is in G) = a m — ^ and hence 



P(K symbol errors in n symbols with s blocks interleaved) 

- a [GG(n,K) + GB(n,K)] + (1 - a) [BG(n,K) + BB(n,K)] 



SECTION IV 
INPUT FOR THE PROGRAM 

XM = No- of bits/symbol 

XN = No. of symbols/block 

XNEST = Largest number of symbol errors to be considered 
(if this field is blank, XNEST = XN) 

XIPER = Number of interleaved symbols (0 or 1 means 1) 

IK = or not equal to (0 means continue with CP, SP, H, SK; 

not equal to means read new parameters) 

CP = P 

SP = p 

H = h 

SK = k 
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SECTION V 

OUTPUT OF THE PROGRAM 

TIMEX = A8 representation of the time read from IBM 7030 Time Clock 
by a STRAP coded routine (one can call his routine or omit 
it altogether) 

CP, SP, H, SK as above 

ALPHA = P (random bit is in G) 

M (or XM) as above 

N (or XN) as above 

NS m No. of recursion terms to be attempted 

WPI - P (symbol error) 

WY155 = P (no errors when s = 1) 

WPMU2 = mean number of errors in a block 

FCMEAN = mean number of errors given an error occurred (when s = 1) 

For J, L, A, B see [2,3,4] 

IS = number of errors 

P = P (IS symbol errors in a block) 

R = P (IS symbol errors in a block |error occurred) 

Q * P (s IS symbol errors in a block) 

QH « P (> IS symbol errors in a block) 

S = mean number of bits between blocks with > IS symbol errors 

PBAR " contribution to mean number of errors per block made by 

probabilities actually calculated (if one wants the whole 
mean, then NEST * or N; the same applies to VAR, SVAR, 
and CMEAN) 
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VAR - contribution to variance about PBAR made by probabilities 
actually calculated by the recursion 

SVAR = approximate standard deviation 

CMEAN = contribution to mean number of errors given an error 
occurred in the block made by probabilities actually 
calculated by recursion 
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APPENDIX 
PROGRAM 



fe^^idS £5gSjf5£ — £OftftECTiatt tEVEt 03/^tf 



OOOOO 
00U01 
00QP2 

0&GO4 
00005 
00006 



-•TO CALCULATE THE RECURSION PROBABI LI T I tS 
COMMON 777/ TIMEX 
DIMENSION GG(515), GB(515), BG(515)t aB(515) 



404 
1701 



«*£*$ttJM B0B*2$5> T 6*»2»* f " BG6t235J f &IM2Q5* 
4»««fiin» £Q£«£5>* £K*Z3SJ, C0Bt205»* G1&U35J 

REA01701, LK 
FOKMAT (12) 






A€JSfr 



00009 
QUO 10 



.REAP IN THE NO. OF BITS PER SYMBOL, THE NO. OH SYMBOLS . 
HO* MANY SYMBOLS TO GO TttftU» Hflk HSfttV 5YH&GLS 

7oi*«**Htn XML** x*e$ 

T (13, IX, 16, 2X 

C. .READ IN CPt SP, H, SK 
604 KEAD 702t IK, CP, SP, H, SK 
702 FOKMAT (II, E17.10, 3(E18.1D) 




iftTERL&AVfc 



500S 

00014 
00015 
00016 


NEST = XNfcST(K) 
IPEKO = XIPfcr<<K) 
IF(NfcST.EQ.O) NEST = N 


--. 



0Q€ 



£ALU UM 

FJtINT 300, TIMfX 




00020 300 FORMAT! 1H1, Ad, 21X, 77H*CALCULA TIUN OF GILBERT CHANNEL ERROR PROBAB 

ULlTIfcS USING RECURSION FORMULAS*/////) 
00021 MUD = M*N 



FLUATtttliO) 

= MUD - 1 
::=— T* - H 




00025 
0CO26 
Q0Q27 



SKCUHP 
CQ = 
SQ = 



1* - SK 

- CP 

- SP 



00031 
00032 
00033 



jGAMHA = L./tCP*5P> 

ALPHA * $P*GArtMA 
„BETA = CP*GA3MA ^. z=f~.-i=£zzr.rE===£=: 

PI = (SP*SKCOMP ♦ CP*HCOMP)*GA*MA 

PRINT 301, CP, SP, H, SK, ALPHA 
301 FORMAT ( 5X, 27HCHANNfcL PARAMETERS 

tlLii ^filW *H w * f F8*^, 5X, 5H K = ptift^lUtSJU &BALPHA 

PRfHT 3Qit ** H* HSi PI* !P€RO 

FORMAT ISX* 2*HE0Q£ PA«A«£TiBS * = 

112X, 35HN0. OF RECURSION TERMS REQUESTfcO 

2 Eli. 7, 75X, 

3 $,/) 



CP = ,E11.6, 5X, 5HSP « 



, 15/ 7X, $P1 = $, 
14, S BLOCKS INTfcRLEAVED 
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0CC3b 
00037 
00036 

00039 
00040 
00041 


wN = 
WT ■ 
HR = 
WbOKT 
wJ = 
WL = 


00042 

00C43 
00U44 
00u45 

30046 


WA * 
1 
WA = 
WH = 
*IX2 = 
WPI 


00047 
00048 


WPMU2 
KJ « 



SP - CQ 

i>K*LQ + H*SQ 

WT*#2 + 4«*H*SK*WN 

= SORT(WR) 

IwT + WS0KT)/2. 

(WT - hSOKTJ/2. 



*(CP*SK*HCOMP 
SP*SKCUM_P ) 



+ SP*H*SKCDHP } / 



WJ + WN 

(CP*HtOMP 4 

WA '- _ lT~ 
( WA*(l.-WJ**M)/( l.-WJJ) - (wB*l l.-WL**M)/t l.-WL J J 
L_ PJ-. * _ wx2 

FLOATtMUO) * ALUUlO(WJ) 



0U049 
00050 
00051 

00052 

00053 
00054 



HL * 

IFtRJ. 

IF*RJ. 

(HRL. 

(FlKL. 

WY155 



fLUAT(MUO) 



LT. 
GE. 
LT. 
GE, 



■50.) 
-50.) 
-50.) 
-50. ) 



ALU&lOtWL) 

= 0. 

k WJ**MUOt^. 

» 0."~" 

= WL**MUD 



00055 
O0056 
£0057- 



l >— f i/ii . - wjn-w d*( i.-o)/( i.-wl n*pi 

1. - WY155) 



= 1. - ( (rtA*( 
f CHfcAW = - WPMU2 / t 

b06f Wt»i ? RTttTf WPMU2 f FCrtfcAN f WJ, WL* WA, Wti 
(IX, 20HPRGB 
,fcl2.6, 1H 
, tl2.6t S*$ / 
HL _= y 3X t _F9 



PR (NT 
JK>6_FGRMAT 
1K0RS = 

2AN = i 
3J5X_, 4 



SYMSflL ERROR = ,£1 2,6« 1H*« 2X* I7HPR0B KO ER 
2X* ah MtAN = ,hl2.6 t IH*, 2X, tCONOCTlUMAL ME 

1JX, 4HJ = , 3X, E16.10, IH*, 
6 t 1H*» 6Xt 4HA = » 3Xt f-9 .6, lH*» 20X f 4HB = 



00055^ 
0C059 



O0C6& 
00061 
J)0062_ 
'00063 
00064 
0CQ65 
00066 
00067 
00068 



4*F9.6i 
"GOGiir 
GlGtl) 

GOBt 1J 
Glb( 1 J 

aobi ij 

eifc(i) 

BOGll) 

_BJG(1) 

IF(M.E 

DO 50 

Jl = J 



iH*///) 

= CG*SK 

^CO^ GQG 

= CP*H 

= CP - GOB 

= SQ*H 

= SO - 




BOB 

BOG 

80 



tl) 
(1) 



0006* 
OO0T& 
00071 



00072 
00073 
00074 



TA 

G16(J)_ 

TB 

GOBI JJ 
GIBt J) 



= SP*SK 

= SP - 
g.n GO TU 

J = 2tM 
-J_ _ 

= G0B(J1)*SP + G0G(J1)*CQ 

* TA*SK 

= TA - GOGtJ) + G1G(JU*CQ + GlSClii^SP 




= G0BU1 )*S0 

= TB*H 

= TB - GOB t J J 



+ GOG( Jl J* CP 
+ G1BU1) *SU + GIGtJl J*CP 



30075 TX "^ B0B(J1J*SQ 

00076 BOBU) * TC*H 

0CO77 B1B U) = TC^ BOB (J) 

00078 TO = B0GU1)*C0 

00079 BUGU) = TO*SK 

00080 50 BIG(J) = TO - BO G! J I 

00081 8 PRINT 303, GOG(M), G1G(H), GOBtH)* GlB(fl) 



♦ B0G(J1)*CP 

♦J5 1 B <J1)*SQ +_ttiAiJU *CP 
+ BOBtJlJ*SP 

+ B1G( Jl)*CQ + BIBtJl J*SP 
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QQOQZ 



303 fORHAr f£'X f 6hG0£ ~ r £12*5* 5XV £H£I£ 



E12.5t 5X, 6HB0B 



000B3 PRINT 304, BOG(M), B1G(M), BOB(M), B1B(M) 

000B4 304 FORMAT <5X, 6HB0G = , E12.b, 5X, 6HB1G ' 

!E12.5t 5X t 6HB1B » , E12.5//) 

Q0085 fcPIl = ALiWA'CGlSfM* * GlGiM) * dElA*<8l8*K* + 5l^(Mn 

OOOS& _- 3»iU&T 555* SPII 

058ST 5£5 f^^MA? iiXj 3Sh**KO& £***OOW SYK8UL IM ERROR » , 



000BB 
000B9 
00090 



155991 

00092 
00093 



CALL TIME 
GG( 1) = GOG(M) 
GG(2) = G1G(M) 
- GAM HI 



" S8f 1J 

GBI2J 
BGMI 



= Gl&ttU 
» 8 0GIMI 




00094 
00095 
00096 
00097 

00099 
00100 
00101 
00102 
GG1G3 



BG(2) = B1G(M) 
BB< 1) = BOB(M) 

BB(2) « B1B(M) 

IF (IPERD .*£• Of 00 10 9001 

Gfl&i = GOG t Ml 

GiGI * GIGJM) 

GOBI = GOB(M) 

G1BI = G1B(M) 

BOG I * BOG(M) 

BiGI - 




00 iv** 
00L&5 



SOtil - Btffi 
B1BI = SS 
SO TO 9 




00107 9001 GXG = l>GU> * GG(2I 
0G1Q8 GXB * GAUI + G&I2I 

Q01Q9 BXG = 6GI1) » BGI2) 



4*0 U0 
00111 



axe * BB( 1) + BBtZl 
TAU ~ GJ£G~i axa - I. 



00113 
00114 
00115 



DEN = 1. - 
GT ■ 1« - 

BT = U - 



TAU 
BXB 
GXG 



00117 



OOUZ 
on 133 
8&L24 

00125 
00126 
00127 



*>€*& = FtftATUPtRflt tt^ zZ__ry 



-^J £ - SS»i 



00119 PRINT 9003, SIG, TAU 

00120 9003 FORMAT ( 8X , SSIGMA = $, E20.12 f 10X» STAU 
00121 AG = tGGtl) *GT + GXB ♦ BG(D) / DfcN 



$, E20.12/* 



• iG^Cl* 



GOGI 
BOGI 

AG ; 



= AG + 

» AB + 

(GG(2) *GT 



CG * G 
CB * G 
+ GXB * BG(2) ) 



/ DbN 



f$Sf^» - AG I I TAU 
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00130 
00131 
00132 
GC133 



CB = *BG(2) - ABJ / TAU 

G1GI = AG + CG * G 

61GI = AB + CB * G 

AG = (GB( 1) *GT + GXtt 



BB( 1)) / OEN 



00134 Cti = (GBUJ 


- AG1 / TAU 


00135 AS = (dXO * 


U3il> 7* $$££) * BTl / Q&i 


00136 C^ = <BB<1J 


— AB J / TAU 


0C137 GCBl = AG 

0013B B0B1 = AB 
00139 AG = (GB(2> 


+ CG * G 
+ CB * G 
*GT + GXB * Bb(<i> ) / OEN 


001*0 CG. = (GB(2) 


- AG J / TAU 


0C141 AJJ = < BXG. * 


Gb*2) * bb(2) * BH / DEN 


00142 CS"r = „"ilB(21_ 





00143 
00144 
00145 



G1B1 = AG + Co * G 
B1B1 = AB + CB * G 
PRINT 900b t G0Gl» GlGlt GOBlt GlBlt BOGlt B1G1, 



BOBIt B1BI 



00146 _ - 9D0* SJflMAl 4/sJ(, 7HG0GI | ,£12^5* £X^ 7HGTGX » 

1 €12^* 4Xf 7H01B1 * * E12*5/ 5X f ThbOgI * 

2 £12-5* 4Xt 7HdOBi ^ « ei£.>» 4X* 7HBIB1 * 



fci2»&f 4X f ?H&i(il j> |g 
£12.5//.! 



00147 9002 "lFtNEST.EQ.l) GO TO 90 

C .CALCULATE THE RECUKS1UN PROBABILITIES 

00146 DO 51 Ll=2tN 



00143 
00150 
00151 
00152 

00153 
00154 



LZ = MINCiLltNESTJ 

LY * CI 

IF4L1 *0T. NEST* LY 



IF1L1 .GT. 

Nl = LI + 
GG(N1) 




NEST) GU TO 97 



GG(L1>*G1GI 



00155 
00156 
00157 

0015b 
00159 
00160 
00161 



G&1N1* 
BG<N1) 

BBtNlJ 



+ GB(LU*B1G( 

GG(L1J*G1BI ♦ £B&1J*SIBI 

6(i<Ll)*GlGl | £fcf*-IJ*£l£T 

BG<L1>*G1BI + fHH*.J**$i$U 



97 CONTINUE 

DU 52 L2=2iLY 
N2 = LY - L2 



+ 2 



N3 = N2 




00162 
00163 
00164 



YAH1 *0G**2J*£0GI 
G3*N2J^GB(NZJ*B0Bi 
£&IHZI = YAtfl 



+GB(N3T*dlBI 



-*Sfr»2i *S$&1 






00165 
00166 
00167 



52 



YAH2 =BG(N2)*G0GI 
BB(N2)=bB(N2)*bOBI 
BG(N2) = YAH2 



+BG(N3>*G1GI 
+BB(N3>*BlbI 



+BB(N2)*bOGI 
+bG(N2)*G0BI 



+BB(N3)*B1GI 
+BG(N3)*G1BI 



001 $S 
00T69 
0C17C 
0G171 
00172 
0C173 



TAH3 

E^^Stl? 

GGi 1} 

YAH4 
BB( n 
BG(l) 



$£t4)*fc0Gl 
GBtlt *BQBi 

YAH 3 
BG( I! 
BttU) 

YAH4 



+ bb{U*ttOGI 
+ Gutl) *GQBI 



*bOGL 
♦BOB! 



btHl) 



♦ BOGI 

♦ GOBI 



00175 
0C176 



?0 



EPP££ 
CONTINUE 




PRINT 
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00177 
00178 
00179 

o&iao 

00181 



503 FURMATi //////////// //50X, 33H( SYMBUL UATA STARTS QN NEXT PACE)) 

IF(NbST.GT.35) PRINT 502 
502 FORMAT ( 1H1) 



00182 
00183 
0Q184 
00185 
5GI86 
00187 



00188 



PH1HT 500, TIMEX 
500 FDRHATrlX»AB/ U* SSYMBUL MUMBtRilOXi ISHftSCURSlG* PftOBS, 3X, 
I^ITHCQNOITIONAL PROBS> ,15X> 13HCUM RECURSION, 12X* .«• r; CUHS, 

BETA *(Bb(l) 4 BB(l) ) 



2 9X, SMEAN •TMEENS) 
P = ALPHA*(GG(1 ) + GB(1) ) 

Q = P 

-fft F~^lj^ P 

51 * * FWO / Rl 
PM1NT 507* Ft 0* Rl* SI 
507 FORMAT (12X,$0*$t 8X, E16.9, 1H*, 34X, E18.12, lH* f 
1 E10.4, $*$, 9X, E10.4, %*%) 
PBAR = 0. 




9X, 



t>0189 


PH- ™3 i~t%$ 


„ 


00191 
00192 
00193 
00194 


P = ALPHA*(GG(i) ♦ GB(I)) + 8ETA *(BG(I) + B6(I>) 

Q = P + 
OH = 1. - Q 




00195 


ft - F / Rl 


— , . -t^t: 


00196 
O019T 


S = a * FHUO / QH 
FIS » FL0AT1IS) 


r— 



00198 FIST= FIS**2 

00199 PBAR = PBAR ^ 

0200 PVAR = PVAR ^ 

QG2G1 ; . PRJNT-501 ,^I:S* -~i 

W20Z 53 CONTINUE 

S02&3 501 FORMAT J BX* 1 5, 

1 E10.4, lh* f 9X, 

00204 606 CONTINUE 

00205 VAR = PVAR - 



00206 
$$20? 
OT20S 
00209 
00210 



P 
Rt 



► FIS 
* FIST 
Qt OH, 



tH* f 2( 9X, El5.8* IM*), 9X t fcl&,U, 
E10.4, $*$) 

PBAR*PBAR 



305 



5 VAR = SQftTfVAfH 
CH£±H - PBAR / Al 

ALL TI«£ 

PRINT 305, P8AR, VAR t SVARt CKEAN. TIMEX 
FORMAT (// 25X, SPARTIAL MEAN 
1 25X f SAPPROXIMATE VARIANCE 



E13. 
E13. 



6/ 
6/ 



tPARTIAL CGKOITIUftAL «£AN ■ *» £13*6// IX, %9/fl 



00212 
00213 
00214 



2SXi 

zsx, 



*TI**UE 

IF(IK .EQ. 
GO TO 404 
END 



0) GO TO 604 
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